Phase diagram of the anisotropic Kondo chain 
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We establish the phase diagram of the one-dimensional anisotropic Kondo lattice model at T — 
using a generalized two-dimensional classical Coulomb gas description. We analyze the problem by 
means of a renormalization group (RG) treatment. We find that the phase diagram contains regions 
of paramagnetism, partial and full ferromagnetic order. 
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The question of the behavior of localized magnetic mo- 
ments in metals bears on a variety of important materi- 
als, from heavy fermion systems to manganites. Central 
to most theoretical studies is the Kondo lattice model 
(KLM) with both antiferromagnetic (AFM) and ferro- 
magnetic (FM) couplings. The former case results from 
superexchange and is associated with heavy fermion be- 
havior in rare-earth and transition metal intermetallic 
systems |Q. On the other hand, direct exchange leads 
to intratomic FM interactions of the Hund's rule type, 
also leading to a KLM description of systems such as the 
manganites ||. In both cases, there is strong interest 
in the determination of the phase diagram and, more im- 
portantly, in the nature of the quantum phase transitions 
that separate the various phases at T = 

Due to the occurrence of non-perturbative Kondo cor- 
relations, in addition to critical order parameter fluc- 
tuations, quantum phase transitions in the AFM KLM 
are still a hotly debated issue [Q] as compared to their 
counterparts in other metallic magnetic systems H (see 
H for a recent attempt). In this letter, we present a 
renormalization group (RG) treatment of the KLM that, 
although confined to one spatial dimension, has the ad- 
vantage of being able to incorporate both types of low- 
energy processes on an equal footing. The RG analysis 
has been very fruitful in the single-impurity case and 
is the adequate tool to address the issue of competing 
ground states. Our treatment generalizes to the lattice 
case the mapping of the single-impurity problem into 
a classical Coulomb gas J?|]. This enables us to deci- 
mate both the conduction electrons and the spins simul- 
taneously. We establish the phase diagram (Fig.|l|) for 
both signs of the coupling constant in a unified fash- 
ion. Our results may prove directly useful for quasi-one- 
dimensional organic compounds with localized moments, 
such as (DMET) 2 FeBr 4 §]. 

The anisotropic KLM chain is described by: 

j'.cr a—x,y 

where c J(T annihilates a conduction electron in site j 
with spin projection a, Sj is a localized spin \ opera- 
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j, the conduction electron 



spin density. At long wavelengths and low energies one 
can linearize the dispersion around the Fermi points ±kp 
(kpa = f n c , n c is the conduction electron number den- 
sity) , and take the continuum limit of the fermionic op- 
erators. Following Ref. || we use bosonization identi- 
ties and neglect the backscattering terms, which are ir- 
relevant away from commensurability. In this limit the 
bosonic charge fields decouple giving rise to gapless col- 
lective modes. On the other hand, the bosonic spin field 
is coupled to the local spins, leading to the Hamiltonian: 

H = T / dx [dx<t>c ' s (a;)]2 + [dx9c ' s (a;)]2 + (1) 
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where Vp = 2t sin kpa is the Fermi velocity and the 
bosonic fields are defined as in Ref. ||. 

As in the Kondo problem, it is convenient to rescale the 
Hamiltonian by the Fermi velocity, introducing the di- 
mensionless coupling constants J z ,± — ^ . We follow 
an approach analogous to the Anderson- Yuval-Haman 
mapping of the Kondo impurity problem onto a classi- 
cal Coulomb gas (CG) j?], [nj. This is achieved by going 
to a path integral formulation in the coherent state basis 
of the bosonic fields and the S z basis of the local mo- 
ments. The z-part is unmodified whereas the transverse 
terms generate spin flips along the Euclidean time direc- 
tion leading to the partition function: 
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Here, y — ^jp, N is the number of flips for a certain spin 
configuration {S z }, and the Euclidean action is: 

S = S + J z ^^2 J dTd x cjy s (x,T)S z (x,T)+ (3) 

-iV^J^ i v (j) 1 U) fa (J) ~ 0s U) 1 V)] > 
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where So is the free Gaussian bosonic action in both vari- 
ables (f> s and 6* s |0|, J2j ls a sum over kink (spin flip) co- 
ordinates and q — S z (r + 5t) — S z (t) = ±1, is called the 
"magnetic charge" [jll| . Because of the cosine in the J± 
term of |l]), for each spin flip the (j> s field comes in with 
two different signs We denote them by v(j) = ±1 

(v(J)q(j) is called the "electric charge"). Thus, each parti- 
cle corresponds to a spin flip and has both a magnetic and 
an electric charge, which are related to the original term 
that produced the flip. The fugacity of the particles is 
y. For instance, a spin flip produced by the right moving 
fermions corresponds to a particle with electric and mag- 
netic charges with the same sign, whereas one produced 
by the left moving fermions gives rise to a particle with 
opposite signs on its charges. There are two restrictions 
on the charge configurations for each space coordinate: 
(i) the magnetic charge q must alternate along the time 
direction (because its origin is a spin- 1/2 flip) and (ii) 
the total charges must be zero 9 = Yl Q v — ® (because 
of periodic boundary conditions in the time direction). 
We note that these conditions are more stringent than 
in usual ID bosonic field theories [|llj. We therefore call 
them strong neutrality conditions. 

The final step consists of tracing out the bosonic fields 
in (H|) in order to obtain an effective action for the spins 
and kinks. When this is done, both short and long range 
interactions are generated. The latter are universal but 
the short range ones depend on the cutoff procedure||l2[. 
These short range terms are essentially the same found 
in reference p^, Q by means of a modified bosonization 
approach. We will focus on the universal long range part 
of the action. Upon integrating by parts in imaginary 
time, spin time derivatives become kink variables. We 
rewrite all long range terms in the form of a generalized 
CG action in two-dimensional Euclidean space for the 
kinks |fjl|: 



S eff = ATlny+fE, 
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where is the Euclidean distance between 2 particles 
and ifiij the angle. In (||), k = 1 — J z /^, 9 = 1, m (j) = 
q (j) is the magnetic charge and e (j) = v (j)q (j) is the 
electric one. The coefficient of the term in ipij is usually 
an integer (the conformal spin, |ll[) and the ambiguity of 
2wn in the angle is then irrelevant. However, in this case, 
k can assume non-integer values. The theory remains 
well defined nevertheless, due to the strong neutrality 
condition, which leads to a cancellation of the Riemann 
surface index. 

In order to investigate the physics of the action M), 
we employ a renormalization group (RG) procedurep5[. 
The most interesting situation is the dense limit, where 
the distance between impurities is of the order of the 
smallest bosonic wavelength available in the system. 



Even though we begin with unitary charges, higher 
charges are generated by renormalization, which come 
from the fusion of elementary kinks ITS]. They correspond 
to new action terms with spin flip pairs and four fermion 
operators: 
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where 8 is a distance of order a. These operators do not 
appear in the original Hamiltonian but are generated by 
the RG procedure. Notice that they are associated with 
exchange processes generated by electron-electron inter- 
actions. It is natural, from this view point, to think of 
the localized spins are generating interactions among the 
electrons. The O p h term flips two nearby spins simulta- 
neously and its action generates a particle with charges 
(m, e) = (±2,0), whereas O pp creates (0,±2) charges. 
Thus, G and G are the fugacities of these charge 2 par- 
ticles. Higher charges are irrelevant. The RG equations 
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where e = sin (2ttk) /(2ttk), with initial conditions: 
g (0) = l,y (0) = Jl/2, and G (0) = G (0) = 0. 

The Coulomb coupling g starts at 1 for non-interacting 
conduction electrons. However, the same RG equations 
apply to the case of conduction electrons with an SU(2) 
non-invariant forward scattering interaction. In this case, 
the initial value of g is the corresponding Luttinger liquid 
parameter|Q. We will not consider this case here, but 
its phase diagram is analogous to the one below. 

By considering the solutions of Eqs. (0), we can trace 
three distinct regions characterized by different fugacity 
flows. In Fig. H these regions are plotted as a function 
of the original Kondo coupling J z and the band filling 
n c . Since the RG equations depend only on \k\, those 
regions are mirror reflections on the n — line. The 
full and dashed lines trace out borders between different 
phases, whereas the dotted line, which is embedded in 
region 3, is the "Toulouse line" of Ref. §. 

In region 1 (k 2 > 3), single spin flip processes are 
irrelevant (y — > 0) just like the FM phase of the sin- 
gle impurity Kondo problemM]. Besides, this phase has 
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Figure 1: Ground state phase diagram of the KLM as a 
function of the Kondo coupling J z and the band filling n c . 
(1) corresponds to the ordered phase of the spin array with 
(S z ) = 1/2, (2) is an ordered phase with (S z ) < 1/2, and (3) 
is a paramagnetic phase, (S z ) = 0. 



one of the higher charges G flowing to strong coupling 
(see ([}])). In contrast, both single and double spin flip 
processes are relevant in regions 2 and 3. What distin- 
guishes them is the fact that in region 2 the flow of y is 
slower and G > G. The flow in the dashed line between 
2 and 3 can be solved analytically: y(£) = y(0)e e and 
G(£) = G(£) = 7n/ 2 (0)(e 2£ - l)/2. This defines a charac- 

1. 



J± 



where y(a/ao) 



teristic length a ~ 2aoln (2/ 
In this case, there is a precise balance between the electric 
and magnetic charges, which prevents them from being 
screened. As a consequence, g = 1 and does not renor- 
malize, though the ground state is a plasma. All cor- 
relations fall in a power law fashion implying a gapless 
system. On the other hand, outside the line \k\ = 1, the 
interactions are screened (g — > or 00). A particularly 
simple case of this kind of flow occurs in the "Toulouse 
line" (k = 0), where g — > and all fugacities grow. 

Although the RG flows are clear, their physical inter- 
pretation is less straightforward. Since we used Abelian 
bosonization, we treated in different ways the z and the 
transverse components of the spins. Therefore, while 
short range transverse spin correlations are generated by 
fusion of elementary particles (Eqs. (||) and the cor- 
responding z correlations appear only through their an- 
nihilation and no fugacity is associated with this pro- 
cess. Nevertheless, we can make progress by writing 
down the operators describing this annihilation. After 
point-splitting the fermionic part, we get: 



O z ~ 2 (G 2 - G 2 ) S" (x + 6) S z (x) - 
,,2 



which are the counterparts of the transverse terms com- 
ing from the fusion of particles. This enables us to deter- 
mine the magnetic phase diagram assigning an effective 
spin Hamiltonian to some special cases. Taking 6 to be 
the lattice spacing at the final RG scale in the O pp , O p h 
and O z definitions, we find an effective Hamiltonian 
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It reproduces the same CG that we studied above. In 
region 1, this reduces to the FM Heisenberg model in 
its ordered phase (G - 1, (S z ) = 1/2). The effective 
Hamiltonian for the n — line is also independent of the 
bosonic field. It is an AFM XYZ model in an external 
field. In this case, the effective spin Hamiltonian is or- 
dered in the XY plane, G ~ G ~ y ~ 1. Nevertheless, 
this does not imply any order of the original spins. A 
unitary transformation connects the spins of Eq. (D to 
the ones of the original model (0) ensuring that the lat- 
ter are disordered even if the former are ordered (see the 
discussion after Eq. ( |Io|) and Ref. iQ). Another situation 
that can be insightful is the |«| = 1 line. In this case, we 
cannot write an effective model for the spins independent 
of the bosonic field. However, due to the symmetric flow 
of G and G, the z-term vanishes and the order parameter 
{S x,y,z ) is still zero. With these assignments in mind, we 
propose that the entire region 3 is a paramagnetic phase 
with short-range AFM fluctuations. There is no simple 
effective Hamiltonian within region 2, but the disorder- 
ing term, proportional to y, starts to grow more slowly 
and the short range z correlations turn from antiferro- 
to ferromagnetic. We thus find that this is an ordered 
phase with unsaturated magnetization of the spins. 

The picture that emerges from these results is that 
there are two continuous phase transitions in the KLM. 
The first transition from region 1 to region 2 in Figjl], 
reminiscent of the Berezinskii-Kosterlitz-Thouless tran- 
sition of the single impurity Kondo model 0, separates 
regions of relevance and irrelevance of the single flip pro- 
cess. The effective model (Eq. ||) in region 1 has FM or- 
der, with full saturation of the localized spins. A regime 
with FM order is beyond the present bosonization treat- 
ment, since the spin polarization of the conduction elec- 
trons must be incorporated. However, the RG flow is still 
able to indicate its existence through the irrelevance of 
single spin flips. In a highly anisotropic model, this leads 
to the ordering of the localized spin array so that the 
electrons can gain kinetic energy (resembling the dou- 
ble exchange mechanism). In Refs. (18, [l9fl, the authors 
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showed that in the isotropic case this is indeed what hap- 
pens. Within this scenario, the total spin per site (elec- 
trons+spins) would be S* ot = (S z ) — n c /2 = (1 — n c )/2 
in the AFM case and Sf ot = (1 + n c )/2 in the FM one. 
Within the region where the couplings are relevant, there 
is another continuous phase transition from region 2 to 
region 3 in Fig.[j], similar to the transition of the Ising 
model in a transverse field [j[7|, that separates a param- 
agnetic phase (region 3 of Fig. |l|) from a region with un- 
saturated magnetization of localized spins, which grows 
continuously until the border of the first transition (re- 
gion 2 of Fig. Q). This interpretation is consistent with 
the numerical studies of both the isotropic FM KLM of 
Dagotto et al. || and the isotropic AFM KLM of Tsunet- 
suguei al. |[xgfl . The methods used here cannot describe 
the region of phase separation found in the numerical 
simulations for the FM KLM, because in that case the 
magnetic energies are of the order of the electron band- 
width. In this limit the bosonization scheme is not ap- 
plicable. 

We now make contact with previous treatments of the 
KLM with Abelian bosonization. In Refs. §, ||, Q, a 
family of unitary transformations, 

is used to define new fields that mix spin and boson de- 
grees of freedom. The charges of the CG we obtained 
arise from vortices of these mixed fields. As can be read- 
ily checked, the integration by parts that we performed 
in the effective action for spins and kinks is equivalent 
to this unitary transformation. Therefore, our effective 
Hamiltonians in the different regions of RG flow (Eq. |) 
should be understood in this rotated basis. This is of spe- 
cial importance in the analysis of the k — line, where 
(9 S ) = 0. Thus, even if the spins acquire order in the 
XY plane, they still remain disordered in the original ba- 
sis. Zachar et al. || argued that along the k — line 
the system has a spin gap. Since the effective Hamil- 
tonian (||) in this line is in a gapped phase, our results 
are consistent with this conclusion. However, this is at 
variance with the available numerical evidence|jl^, 
for the isotropic KLM. This discrepancy raises the ques- 
tion about whether the anisotropic model can capture 
the physics of the isotropic one. In addition, the sub- 
sequent work of Zachar proposes additional phases 
away from k = 0, which may be related to our regions 1 
and 2. Honner and Gulacsi[jl3|, have also proposed a 
phase diagram for the isotropic ID KLM. They predict a 
paramagnetic phase for ferromagnetic coupling. This is 
in disagreement with our results and the work of Dagotto 
et a/|Q. A full discussion of their methods and results in 
contrast to ours will be published elsewhere. 

In conclusion, we have established the zero temper- 
ature phase diagram of the anisotropic ID KLM with 
ferromagnetic and antiferromagnetic coupling. We have 



found three different phases: a paramagnetic phase where 
the Kondo effect dominates; a fully polarized magnetic 
phase where the "double exchange" correlations drive the 
system towards order; and a partially polarized phase 
where Kondo effect and magnetic correlations compete 
directly to generate partial polarization. The two quan- 
tum phase transitions have continuous nature, closely 
related to the Berezinskii-Kosterlitz-Thouless transition 
of the single impurity Kondo problem and to the Ising 
model in a transverse field. Although we have worked in 
ID, many of the effects discussed here are generic and 
also occur in higher dimensions. In spite of the fact that 
we have used Abelian bosonization and worked on the 
anisotropic model, our findings are in agreement with the 
numerical simulations in the SU(2) KLM [HQ. It would 
be interesting to have the numerical work extended to the 
anisotropic model as a further test of our results. Finally, 
we hope our results will be a stimulus for the study of 
the phase diagram of quasi-one-dimensional systems with 
localized moments such as (DMET)2FeBr4 ||, 
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